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^ ■ Abstract 



We present a calculation of two photon radiation in W and Z boson pro- 
5^ \ duction in hadronic collisions, based on the complete matrix elements for the 

processes qq' — > £^i"yj and qq — > i^i~'yj, including finite charged lepton 
masses. In order to achieve stable numerical results over the full phase space, 
multiconfiguration Monte Carlo techniques are used to map the peaks in the 
differential cross section. Numerical results are presented for the Fermilab 
Tevatron. 
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I. INTRODUCTION 

The Standard Model of electroweak interactions (SM) so far has met all experimental 
challenges and is now tested at the 0.1% level [|l|. However, there is little direct experimental 
information on the mechanism which generates the masses of the weak gauge bosons. In the 
SM, spontaneous symmetry breaking is responsible for mass generation. The existence of 
a Higgs boson is a direct consequence of this mechanism. At present the negative result of 
direct searches performed at LEP2 imposes a lower bound of Mh > 98.8 GeV on the Higgs 
boson mass. Indirect information on the mass of the Higgs boson can be extracted from the 
Mh dependence of radiative corrections to the W boson mass, Mw, and the effective weak 
mixing angle, sin^ O^Fj . Assuming the SM to be valid, a global fit to all available electroweak 
precision data yields a (one-sided) 95% confidence level (CL) upper limit on M^ of about 
260 GeV [IHI. 

Future more precise measurements of Mw and the top quark mass, m^p, will lead to 
more accurate information on the Higgs boson mass @-§[. Currently, the W boson mass 
is known to ±42 MeV [^ from direct measurements. The uncertainties of the individual 
experiments contributing to this value are between about 80 MeV and 110 MeV 0J^. The 
present uncertainty of the top quark mass from direct measurements is ±5.1 GeV |^. With 
a precision of 30 MeV (10 MeV) for the W mass, and 2 GeV for the top quark mass, 
Mh can be predicted from a global analysis with an uncertainty of about 30% (15%) [^If. 
Comparison of these indirect constraints on Mh with the results from direct Higgs boson 
searches at LEP2, the Tevatron collider, and the Large Hadron Collider (LHC) will be 
an important test of the SM. They will also provide restrictions on the parameters of the 
Minimal Supersymmetric extension of the Standard Model (MSSM) [|T0[] . 

A significant improvement in the W mass uncertainty is expected in the near future 
from measurements at LEP2 |Tl| and the Fermilab Tevatron pp collider [Q. The ultimate 



precision expected for Mw from the combined LEP2 experiments is 30 - 40 MeV [^. At 
the Tevatron, integrated luminosities of order 2 fb^^ are envisioned in the Main Injector 



Era (Run II), and one expects to measure the W mass with a precision of approximately 
40 MeV P per experiment. The prospects for a precise measurement of Mw would further 
improve if a significant upgrade in luminosity beyond the goal of the Main Injector could be 



realized. With recent advances in accelerator technology |T2[ , Tevatron collider luminosities 
of order 10'^^ cm^^ s^^ may become a reality, resulting in integrated luminosities of up to 
10 fb~^ per year. With a total integrated luminosity of 30 fb~^, one can target a precision 
of the W mass of 15 - 20 MeV [^ . A similar or better accuracy may also be reached at the 
LHC [IH]. 



In order to measure the W boson mass with high precision in a hadron collider envi- 
ronment, it is necessary to fully understand and control higher order QCD and electroweak 
(EW) corrections to W production. The determination of the W mass in a hadron col- 
lider environment requires a simultaneous precision measurement of the Z boson mass, Mz, 
and width, Tz- These quantities serve as reference points. When compared to the value 
measured at LEP, they help to accurately determine the energy scale and resolution of the 
electromagnetic calorimeter, and to constrain the muon momentum resolution ^. In order 
to extract M\y from hadron collider data, it is therefore also necessary to understand the 
higher order QCD and EW corrections to Z boson production in hadronic collisions. 

Electroweak radiative corrections have a significant impact on the W and Z boson masses 
and widths extracted from experiment. Recent improved calculations of the 0{a) EW 
corrections to W production |T^, and of the 0{a) QED corrections to Z production in 



hadronic collisions [|T^, have shown that the main effect is caused by final state photon 



radiation. When detector effects are included, 0{a) radiative corrections shift the W mass by 
about —50 MeV in the electron case, and approximately —160 MeV in the muon case |[T6| , p!7| . 
The effect on the Z mass is about a factor two larger than that on M]^ for both electron 
and muon final states. 0{a) photon emission also shifts the width of the W boson extracted 



from the tail of the transverse mass distribution by approximately —70 MeV |[T^. The size 
of the shift in Myy, Mz and the W width introduced by the 0{a) corrections raises the 
question of how strongly (9(a^) corrections affect these quantities. 



In order to reliably calculate the impact of the 0(q;^) corrections to pp —>■ W 
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and pp — * 7*, Z — !► i'^i~ (i = e, //) on the W and Z masses extracted from experiment, a 
full calculation including real and virtual corrections, which is valid over the entire allowed 
phase space, is needed. So far, only partial calculations for the O(a^) real photon corrections, 
pp — ^ W^ -^ i'^z/77 and pp ^ 'j*,Z ^ £+£"77 exist [|l^,Q. In Ref. [|l^] the structure 
function approach for photon radiation is used to perform the calculation and, therefore, 
only final state photon radiation in the leading log approximation is included. The results 
obtained using this approach are reliable only for small opening angles between the photons 
and the charged leptons. The calculation of Ref. pO[ is based on the full set of tree level 



Feynman diagrams contributing to £z/77 and £"'"£"77 production. In addition, to preserve 
gauge invariance when finite W width effects are included, the imaginary part of the WW'j 
and WW'j'j one- loop vertex corrections is taken into account. However, charged leptons are 
assumed to be massless, and thus a finite lepton - photon separation cut has to be imposed 
in order to avoid the coUinear singularities associated with final state radiation. 

The first step towards a calculation of the (9(a^) corrections to W and Z boson produc- 
tion in hadronic collisions thus is to perform a calculation of iwyj and i^i~'-f'-f production 
which 

• is based on the full set of Feynman diagrams contributing at tree level, 

• includes finite lepton mass effects, 

• is gauge invariant when finite W width effects are taken into account, 

• and is valid for arbitrary lepton - photon opening angles. 

In addition, in order to obtain reliable information on the shift in Mw and Mz caused by two 
photon radiation, the numerical calculation should be stable for photon energies as small as 
the tower threshold of the electromagnetic calorimeter of the Tevatron experiments, which 
is of (9(100 MeV). In this paper we present such a calculation. 



While the calculation of the qq -^ £"'"£"77 and qq' -^ ivyy matrix elements is straight- 
forward, the phase space integration presents some challenges, due to the sharp peaks in 
the matrix elements which arise from the soft and coUinear singularities. The coUinear 
singularities associated with final state radiation are regulated by the finite mass of the 
leptons whereas soft and initial state coUinear singularities are rendered finite by transverse 
momentum cuts imposed on the photons. Both soft and coUinear singularities produce 
large contributions to the cross section in small regions of phase space. Standard adaptive 
Monte Carlo integration routines such as VEGAS |2l]| do not yield a numerically stable 
result of the cross section for processes which exhibit a complicated peaking structure in 
the matrix elements. To obtain numerically stable and accurate results for such processes, a 
multi-channel Monte Carlo approach p2|, augmented by the adaptive weight optimization 



procedure described in Ref . p3| , is frequently used. The disadvantage of the multi-channel 



Monte Carlo approach is that the peaks in the matrix elements have to be mapped by hand, 
thus requiring a substantial amount of analytic work which has to be repeated for each new 
process one wishes to analyze. 

Our calculation is based on a similar approach which adds the benefit of largely automiz- 
ing the mapping of the peaks in the matrix elements. The process independent features of 
our approach, and the resulting multiconfiguration Monte Carlo (MCMC) integration pro- 
gram, are briefly described in Sec. II. Full details will be given elsewhere [^]. In Sec. Ill we 
discuss technical detaUs associated with the calculation of the qq — > £+£^77 and qq' -^ ^z/77 
matrix elements and present numerical results for two photon radiation in W and Z events 
at the Tevatron collider {pp collisions at 1.8 TeV). Finally, summary remarks are given in 
Sec. IV. 

II. PHASE SPACE INTEGRATION 

The matrix elements for qq -^ £"'"£"77 and qq' — > £z/77 have many sharply peaked regions 
throughout phase space. In addition to the Breit-Wigner resonances around the W or Z pole 



and a pole at small £"'"£" invariant masses due to photon exchange, there are singularities 
when either photon becomes soft, or collinear with a charged particle. Although these soft 
and collinear singularities are regulated by energy or transverse momentum cuts and fermion 
masses, they result in large contributions to the cross section over relatively small regions 
of phase space and cause difficulties for standard integration techniques. 

Integrating over the parton distributions and final state momenta in general requires 
performing a (3A^/j„a«— 4)+2 dimensional integral over a phase space which may include many 
cuts. Here N final is the number of particles in the final state. If the number of dimensions 
is large, the integral is most easily carried out using Monte Carlo techniques. Monte Carlo 
integration approximates the integral by taking the average of a number of points, A^, selected 
at random, and multiplying by the volume, V , over which one is integrating, 

Jf(x)dxC:^j^J2fix,)xV. (1) 

Provided the function f{x) which is to be integrated is sufficiently fiat, the number of points 
for convergence is independent of the number of dimensions. However, if f{x) is sharply 
peaked convergence may be exponentially slow. 

In order to use Monte Carlo techniques for integrating a sharply peaked function, it is 
necessary to remove the peaks. Peaks which are analytically integrable, and for which the 
integral is invertible, can be smoothed with the appropriate transformation of variables. 
A Breit-Wigner resonance is an excellent example for such a case. The transformation 
y = arctan(x) removes the peak and makes the integrand fiat. Collinear and soft poles often 
require more involved transformations which may not have general analytical solutions. 

Adaptive Monte Carlo programs such as VEGAS are able to fiatten peaks by using 
numeric approximations of the integrand. The result is not as fast, or efficient as analytically 
removing the peaks, however it is more convenient. For most applications this is very 
desirable. The major restriction is that programs such as VEGAS can only remove peaks 
which are in the plane of one of the integration variables. For example, these programs will 
successfully fiatten the peaks in the function 
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f{x,y) = --, (2) 

X y 

but they will not be able to flatten those for 

f{x,y) = — ■ , (3) 

X + y X — y 

unless a change of variables is performed. For processes with relatively few peaks, it is 
usually possible to map each peak to one of the integration variables. For complicated 
process such as qq -^ £"'"£"77 and qq' -^ £h"yy, this is not the case. 

In cases where it is not possible to simultaneously map every peak to an integration 
variable, there are two classes of solutions available. The first is to divide up phase space 
with cuts, such that the peaks in each region can be mapped to the integration variables. 
An adaptive Monte Carlo integration routine is used for each region separately. The results 
from each region are combined to obtain the total cross section. This method is effective for 
processes with relatively simple peaking structure, however as the number of peaks increases 
the technique quickly becomes cumbersome and prone to error. 

The second technique P^J23[] is to choose points in phase space not according to a sin- 



gle distribution, but according to the sum of multiple distributions. Each distribution is 
responsible for a specific set of peaks. The optimal number of points from each distribution 



is chosen using an algorithm which minimizes the Monte Carlo integration error [^. With 
each channel, a different set of poles is analytically removed. The resulting code is very fast 
and efficient, however it requires significant analytic work, which must be repeated for each 
new process. 

In our approach, we have combined the power of multichannel integration with the 
convenience offered by adaptive Monte Carlo integration routines such as VEGAS. The 
result is a general and flexible multiconflguration Monte Carlo program called MCMC which 
can numerically integrate sharply peaked functions in many dimensions with minimal input 
from the user. 

Feynman diagrams offer a convenient mechanism for determining in which dimensions 
peaks may appear. At tree level strong peaks in the cross section are always associated with 



a propagator going on-shell. We have implemented a general phase-space generator based on 
Feynman diagrams. Given a tree-level Feynman diagram, it maps a set of random numbers 
to a point in phase space such that each propagator represents one of the dimensions of 
integration. The operation is invertible so it can also return the set of random numbers 
associated with any point in phase space. The user specifies the momentum flow of the 
contributing Feynman diagrams in a simple include file, together with the masses and widths 
of the Breit-Wigner resonances which appear in each diagram. All other aspects of the phase 
space integration are handled automatically by the program. A more detailed description 
of the approach will be given elsewhere [^ . 



The convenience of MCMC is best illustrated in a simple example. Consider the process 
i^fi^^fj. ~^ e^e~77 for a center of mass energy of ^/s = 100 GeV, where each of the final state 
particles is required to have a transverse momentum px > 10 GeV. The electron mass is 
assumed to be variable. The six Feynman diagrams associated with this process can easily 
be generated with a program such as MadGraph [^ . The diagrams generated by MadGraph 
are shown in Fig. |T[ Each diagram represents a phase space configuration in the MCMC 
code with the appropriate poles mapped to the integration variables. 

These configurations are input to the integration package, which then searches for peaks, 
and determines the optimal number of points to choose from each configuration using an 
algorithm which minimizes the integration error. Table | compares the cross sections ob- 
tained using the traditional single configuration approach with those from MCMC as the 
electron mass is varied from 0.01 GeV to 10 GeV. Notice that for large masses, the matrix 
element is relatively flat and a single configuration accurately integrates the cross section. 
However, as the electron mass decreases, the contribution from the collinear regions becomes 
increasingly important, and the single configuration package is unable to accurately inte- 
grate the function. Not only is the error larger, but even with 5 x 10^ integration points, the 
single configuration integration is giving the wrong result as it samples all of the points from 
the peaks it has mapped to integration variables, completely neglecting the peaks which 
are only sampled by the other configurations. Due to the mass singular terms associated 
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with final state radiation in the colhnear hmit, the cross section scales approximately with 
{\og{ml/s)y for small electron masses, rrie < 1 GeV. This provides a simple check on the 
accuracy of the MCMC result. 

The primary difference between our approach and other multichannel techniques is its 
generality. The mapping of uniformly distributed random numbers to points in phase space 
can be broken down into two steps. First the uniformly distributed random numbers are de- 
formed into non-uniform numbers, with a corresponding Jacobian. Next these non-uniform 
numbers are mapped to four-momenta in phase space with another Jacobian. For each 
step, one can choose to perform an analytic transformation which will be very efficient for 
the particular process being studied, or one can choose a general transformation which is 
not optimized for the specific process, but will work for any process. The approach of 
Ref. ||23| produces highly optimized code for both transformations. In Ref. [^, a general 
procedure for the transformation from uniform space to a deformed space is given, but 
optimized procedures for the transformation to phase space are chosen. The MCMC pro- 
gram provides general algorithms for both transformations, similar to the approach used by 
CompHEP [0,0 to perform the phase space integration. The resulting code is in general 
slightly slower than that resulting from the other two approaches, however we believe its user 
friendliness makes up for this short coming. The advantage of user friendly programs at the 
expense of computer time has already been demonstrated by packages such as MadGraph, 



CompHEP p7| , |29[| and GRACE |2^ which quickly produce non-optimized tree-level matrix 
elements. Indeed the synthesis of the integration package outlined here with automatically 
generated matrix elements will allow the user to concentrate on the physics issues rather 
than numerical integration techniques. While MCMC naturally interfaces with MadGraph 



and HELAS |30|, matrix elements resulting from any other automated or non-automated 



calculation can be used. 



III. tv-i-i AND ^+^-77 PRODUCTION AT THE FERMILAB TEVATRON 

We shall now discuss the calculation of £z/77 and £^£"77 production in hadronic colli- 
sions, together with some phenonienological applications relevant for future W mass mea- 
surements at hadron colliders. To calculate the matrix elements for qq — * £+£^77 and 
qq' -^ i^w-f-f we use MadGraph which automatically generates the SM matrix elements 
in HELAS format. When photon exchange is taken into account, 40 Feynman diagrams 
contribute to £^£^77 production, while there are 21 diagrams for i^U'yy production. Tak- 
ing into account symmetries in the phase space mapping, 20 (12) different configurations 
contribute to qq -^ £+£"77 {qq' -^ iu'j'j). 

In order to maintain electromagnetic gauge invariance for qq' -^ i^wy'y in presence of 
finite W width effects, the W propagator and the WW'j and WW'j'j vertex functions in 
the amplitudes generated by MadGraph have to be modified [^,0. Finite width effects 
are included by resumming the imaginary part of the W vacuum polarization, Uiylq"^). The 
transverse part of Uwiq"^) receives an imaginary contribution 

lmllliq^)=q^^ (4) 

while the imaginary part of the longitudinal piece vanishes. The W propagator is thus given 
by 



D'wiQ) 



q^-M^ + iq^^w 



q^'q" 



9"" - \^a + ^iw) 



M^ 



(5) 



with 



where Tw denotes the W width. A gauge invariant expression for the amplitude is then 
obtained by attaching the final state photons to all charged particle propagators, including 
those in the fermion loops which contribute to Uwiq"^)- As a result, the lowest order WW'y 



"/^Ai r^r^A r'^l^t^p 



and WW'yy vertex functions, Fg and Fg ^^, are modified p0|j31[] to 
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T'^^^ = Tf^{l + ilw). (7) 

T''^^'' = Tf^''{l + i^w). (8) 

The SM parameters used in our numerical calculations are My/ = 80.3 GeV, Tyr = 
2.046 GeV, Mz = 91.19 GeV, Tz = 2.49 GeV, and a(M|) = 1/128. These values are 
consistent with recent measurements at LEP, LEP2, the SLC and the Tevatron |^. We use 
the parton distribution functions set A of Martin- Roberts-Stirling [^ with the factorization 
scale set equal to the parton center of mass energy a/J. All numerical results are obtained for 
pp collisions with a center of mass energy of ^/s = 1.8 TeV. In Run II, the Tevatron collider is 
foreseen to operate at y/s = 2 TeV. For a center of mass energy of 2 TeV, results qualitatively 
similar to those reported here are obtained. Cross sections are about 5% higher than those 
found for y/s = 1.8 TeV. Since the total cross sections for £+z/77 and i'U'yy production are 
equal in pp collisions, we shall not consider the i'wyj channel in the following. 

To simulate the fiducial and kinematic acceptances of detectors, we impose the following 
transverse momentum (pt) and pseudo-rapidity [r]) cuts on electrons and muons: 

electrons muons 

Pj^ie) > 20 GeV p^ifJ-) > 25 GeV 

\v{e)\ < 2.5 |7?(/x)| < 1.0 

l^j, > 20 GeV ^j, > 25 GeV 

Here, ^rp denotes the missing transverse momentum which we identify with the transverse 
momentum of the neutrino in iw-f-f production. The ^j, cut is only applied in iwyj produc- 
tion. The cuts listed above approximately model the acceptance of the CDF detector for 
electrons and muons in Run I. Qualitatively similar numerical results are obtained if cuts 
are used which approximate the phase space region covered by the upgraded CDF detector 



for Run II |33|, or if cuts are used which model the acceptance of the D0 detector |34 



In addition to the lepton cuts listed above, a pseudo-rapidity cut 

1^(7)1 < 3.6, (9) 
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and a transverse momentum cut on the photons are imposed. In order to be able to ac- 
curately determine the shift in the W and Z boson masses induced by photon radiation 
correctly, it is necessary to consider photon transverse momenta as low as the calorimeter 
threshold of the detector, which is about 100 MeV. Subsequently, we therefore require 

Pt(7) > 0.1 GeV (10) 

in all our calculations unless stated otherwise explicitly. The photon transverse momentum 
and pseudo-rapidity cuts are necessary to avoid soft singularities and coUinear divergences 
associated with initial state radiation. 

Since we are mostly interested in photon radiation in W and Z decays, we impose 
additional cuts on the di-lepton invariant mass, 

75 GeV < m{ee) < 105 GeV, (11) 

and the transverse mass of the iu system, 

65 GeV < mr(%) < 100 GeV. (12) 

CDF and D0 utilize similar cuts in their W mass analyses [|T^,0. Events satisfying Eqs. (|TTp 
and (|l^) are called Z -^ i^i~ and W ^ iu events, respectively, in the following. 

To demonstrate that the multiconfiguration Monte Carlo approach we use yields accurate 
results both in the coUinear region as well as for photons emitted at large angles, we show 
in Fig. 1^ the differential cross section versus the separation between the two photons in the 
azimuthal angle- pseudorapidity plane. 



AR,, = ^A02^ + Ar/2^ . (13) 

The strong peak for small AR^^ arises when both photons are emitted by the same charged 
lepton, and the photons are coUinear with the lepton. The peak at AR-y^ ?« 3 in £"'"£"77 
production originates from Feynman diagrams where the photons are radiated off different 
leptons. Since photons do not couple to neutrinos, this peak is absent in lu^yy production. 
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For electrons, the collinear peaks are significantly more pronounced than for muons. The 
difference in the differential cross sections for electrons and muons away from the collinear 
regions is entirely due to the different pt and rapidity cuts imposed on these particles. The 
statistical fluctuations are quite uniform over the full range of AR^.y values considered, 
indicating that the MCMC program distributes the generated events and their weights 



appropriately. Away from the collinear peaks, our calculation^ agrees with that of Ref. ^^ 
to better than 1%. In this region, conventional adaptive Monte Carlo routines such as 
VEGAS are sufficient in order to obtain a numerically stable result. The distributions of 
the separation between the photons and the charged lepton (leptons) are qualitatively very 
similar to the AR^^ spectrum. 

In Tables H and pT| , we list the fraction oi W ^ iu and Z -^ i~^i~ events at the Tevatron 
which contain two photons as a function of the minimum photon transverse momentum. For 
comparison, we also list the event fractions containing one photon. Fractions are obtained 
by normalization with respect to the lowest order cross section within cuts. The results for 
iwy and i'^i^'y production are obtained using the calculation of Ref. PT|. No lepton-photon 
or photon-photon separation cuts are imposed. 

Approximately 3% (1%) of all W ^ eu {W ^ fxu) events, and 14% (5%) of all Z -^ 
e'^e" {Z — > fi^fi~) events, contain two photons with a minimum transverse momentum of 
^^"■(7) = 0.1 GeV. Because of the mass singular logarithms associated with final state 
photon bremsstrahlung in the collinear limit, the fraction oi W ^ eu and Z — > e^e~ events 
with two photons is more than a factor 3 larger than the corresponding fraction oi W ^ fiu 
and Z — i> fi~^fi~ events. In contrast to W events, both leptons can radiate photons in Z 
decays. As a result, the probability oi Z ^ £^i~ events to radiate two photons is more than 
four times that oi W -^ iu events. For increasing p™"'(7), the fraction of W and Z events 



^Parton level FORTRAN programs for pp ^ £ 1/77 and pp -^ i~^£ 77 which include the MCMC 
source code are available upon request from the authors. 
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containing photons drops quickly. 

For small photon transverse momenta, the cross section is completely dominated by final 
state radiation. In this region, the fraction of events containing two photons, P2, can be 
estimated using the simple formula [^ 

^2 = f , (14) 

where Pi is the fraction of events containing one photon. The results obtained using Eq. ([Mf) 
are also listed in Tables || and |T|. For large values of pr(7), the available phase space for final 
state radiation is strongly reduced by the transverse momentum cuts imposed on the leptons, 
and initial state radiation plays an increasingly important role. For p^*"(7) > 3 GeV, 
Eq. (|14D therefore becomes more and more inaccurate. 

The large mass singular terms associated with final state bremsstrahlung result in a 
significant change in the shape of the mT{i^j') and the di-lepton invariant mass distributions. 
This is demonstrated in Fig. ^ Here we do not impose the di-lepton invariant mass cut and 
the iu transverse mass cut of Eqs. ( ]TT| ) and (|T^). In Fig. ^ we show the ratio of the £"'"£"77 



and the lowest order i'^i~ cross section as a function of m{ii). The cross section ratio is 
seen to vary rapidly. The dip at m{ii) = Mz is a direct consequence of the Breit-Wigner 
resonance of the Z boson. Below the Z peak, the cross section ratio rises very sharply and in 
the region 70 GeV < ■m(ee) < 80 GeV, the cross section ratio is of order one in the electron 
case. The dip located at m{ii) = Mz and the substantially enhanced rate of events with 
two photons below the resonance peaks are caused by final state bremsstrahlung in events 
where the £"'"£"77 invariant mass is close to Mz- 

Fig. Pd displays the ratio of the i'^w-f-f and the i'^u cross section as a function of the 
iu transverse mass. Here the dip at mj'{i'^rp) = My/ is due to the Jacobian peak in the lu 
transverse mass distribution. Because of the long tail of the lowest order mTii^rp) distribu- 
tion below Mw and the fact that photons are not radiated by neutrinos, the enhancement in 
the £1/77 to £z/ cross section ratio is less pronounced than that encountered in the i~^i~ case. 
In the region of large transverse masses, mT{i^j') > 100 GeV, the shape of the transverse 
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mass distribution is sensitive to the W width. Fig. ^d shows that two photon radiation 
significantly modifies the shape of the ■mj'{i'^rp) distribution in this region. This will directly 
infiuence the W width extracted by experiment. 

The shape changes in the iu transverse mass and the di-lepton invariant mass distribu- 
tions suggest that two photon radiation may have a non-negligible effect on the measured 
W and Z masses, and also on the W width extracted from the high transverse mass region. 
Since the shape change caused by two photon radiation in the distribution used to extract 
the mass is more pronounced in the Z case, the shift in the Z boson mass is expected to be 
considerably larger than the shift in Mw- For a realistic calculation of how 0(«^) correc- 
tions affect the W and Z resonance parameters, soft and virtual corrections and detector 
resolution effects need to be included. 

We have not taken into account detector resolution effects or realistic lepton and photon 
identification requirements in the calculations presented in this Section. In particular, we 
have assumed that photons and leptons with arbitrary small opening angles can be discrim- 
inated. In practice, the finite resolution of the electromagnetic calorimeter makes it difficult 
to separate electrons and photons for small opening angles between their momentum vectors. 
Electron and photon four-momentum vectors are therefore recombined if their separation 
in the azimuthal angle-pseudorapidity plane is smaller than a critical value [Il6| , |l7| . This 



eliminates the mass singular terms associated with final state photon radiation and thus 
may reduce the fraction of W and Z events with two photons significantly. Since muons 
are identified by hits in the muon chambers, the four momentum vectors of muons and 
photons are not combined for small opening angles. Instead, one frequently requires the 
photon energy to be below a threshold Ec in a cone around the muon. The mass singular 
logarithms thus survive in the muon case. The precise lepton identification requirements 
and their effects on the size of the EW corrections W and Z boson production are detector 
dependent. 
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IV. SUMMARY AND CONCLUSIONS 

The mass of the W boson is one of the fundamental parameters of the SM and a precise 
measurement of M^ is an important objective for current experiments at LEP2 and future 
experiments at the Tevatron. A precise measurement of Mw helps to constrain the Higgs 
boson mass from radiative corrections. It will also provide restrictions on the parameters 
of the MSSM. In order to perform such a measurement at a hadron collider, it is crucial 
to fully control higher order QCD and EW corrections to W production. In a precision 
measurement of M^ in hadronic collisions, a simultaneous determination of the mass of the 
Z boson is required for calibration purposes. A detailed understanding of the QCD and 
electroweak corrections to Z boson production is therefore also necessary. 



Recent calculations ||l^ , p!5[| have shown that the 0{a) electroweak corrections to W and Z 
production have a significant impact on the weak boson masses extracted from experiment. 
The dominant contribution originates from final state photon radiation. The magnitude of 
the shift in Mw and Mz induced by the 0{a) corrections suggests that 0(q;^) corrections 
may have an effect which cannot be ignored in future W mass measurements at the Tevatron. 
In this paper we have presented a calculation of the real 0{a'^) photonic corrections to W 
and Z boson production in hadronic collisions. Our calculation is based on the full set of 
Feynman diagrams contributing to £^£"77 and ivyy production and includes finite lepton 
mass effects. In order to maintain gauge invariance in £1/77 production, the W propagator 
and the WW'j and WW'j'y vertex functions are modified using the prescription given in 
Refs. H^ and ill . 



In order to accurately determine the shift in the W and Z masses caused by photon 
radiation, the numerical calculation should be stable for arbitrarily small or large lepton 
- photon opening angles as well as for photon energies as small as the tower threshold of 
the electromagnetic calorimeter of the Tevatron experiments, which is of (9(100 MeV). Due 
to the collinear and soft singularities present, this poses a challenge. Standard adaptive 
Monte Carlo integration routines such as VEGAS do not yield a stable result for processes 
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with a complicated peaking structure in the matrix elements, such as qq -^ £"'"^"77 and 
qq' —>■ iw-f-f. To obtain numerically stable and accurate results in these cases, multi-channel 
Monte Carlo integration techniques are frequently used. This approach requires that the 
peaks in the matrix elements are analytically mapped. To calculate the cross sections for 
£"'"£"77 and iwy'y production at hadron colliders, we developed a multiconfiguration Monte 
Carlo integration routine called MCMC which is based on a similar approach, adding the 
benefit of largely automizing the mapping of the peaks. MCMC thus can be used to calculate 
other processes with matrix elements exhibiting a complex set of peaks with almost no 
additional effort. The algorithm which is used in MCMC to map out the peaks is based on 
the Feynman diagrams which contribute to the process considered. 

Imposing W and Z boson selection cuts on the final state leptons, we found that a 
significant fraction of weak boson events contains two photons. The probability for Z events 
to radiate two photons is almost a factor five larger than that for W events. For W ^ fiu 
and Z —>■ fi~^fi~ decays, the rate for two photon radiation is about a factor 3 smaller than 
the corresponding rate for decays with electrons in the final state. If the photon pt is less 
than about 3 GeV, the fraction of W and Z events containing two photons can be estimated 
with an accuracy of 20% or better using a simple equation (see Eq. (|T^). 

Two photon radiation was also found to significantly alter the shapes of the Z boson 
resonance curve and the iu transverse mass distribution. The shift in the W and Z masses, 
and in the W width measured from the tail of the transverse mass distribution, caused 
by the O(a^) real photon corrections may thus be non-negligible for future hadron collider 
experiments. For a realistic estimate of how strongly the (9(a^) corrections affect the W 
boson parameters extracted from experiment it is necessary to include the effects of soft 
and virtual corrections, as well as detector resolution effects. The calculation of £"'"£"77 and 
iw-ff production presented in this paper thus only is the first step towards a more complete 
understanding of the O(a^) electroweak corrections to W and Z production in hadronic 
collisions. 
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TABLES 
TABLE L Integrated cross section for the process f^P^ -^ e^e~^^ at y/s = 100 GeV as a 
function of the electron mass for single and rtiulti configuration adaptive Monte Carlo integration. 
A pt > 10 GeV cut is imposed on all final state particles. In all cases 8 x 100, 000 events are 
generated to set the grid, and 5x1 million events for evaluating the integral. 



electron mass 

(GeV) 



Multiconfiguration 
a(fb) 



X 



Single configuration 
a(fb) 



X 



0.01 

0.1 

1.0 

10.0 



88.4 ± 0.8 0.2 

46.34 ± 0.08 0.7 

17.14 ± 0.03 1. 

1.999 ± 0.003 0.3 



48.9 ± 3 0.3 

36.5 ± 2 1.1 

16.5 ± 0.2 0.7 

1.992 ± 0.003 0.2 
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TABLE II. Fraction ofW^^ev and W ^>- fiv events (in percent) containing one or two photons 
with a transverse momentum pt{i) > ^^^(7) at the Tevatron {pp cohisions at y/s = 1.8 TeV). 
Fractions are obtained with respect to the lowest order cross section. The cuts imposed are specified 
in the text. The relative statistical error on the event fractions from the Monte Carlo integration 
is approximately 1%. 



p^*"(7) (GeV) 



W — > ei/7 



W -^ ei/77 



W -^ ei/77 [Eq. (0)] 



0.1 
0.3 

1 

3 
10 



p^*"(7) (GeV) 



23.9 

17.3 
10.4 
4.82 
0.56 



W — > /xz/7 



3.05 
1.56 
0.53 
0.09 
1.3 X 10-3 



W — > /izv77 



2.86 
1.50 
0.54 
0.12 
1.6 X 10-3 



W ^ /.Z.77 [Eq. (|T§] 



0.1 
0.3 

1 

3 
10 



13.5 
9.65 
5.74 
2.65 
0.33 



0.99 

0.48 

0.17 
2.7 X 10-2 
6.3 X 10"^ 



0.91 

0.47 

0.16 
3.5 X 10-2 
5.4 X 10"^ 
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TABLE III. Fraction of Z ^> e~^e~ and Z -^ /^^A*" events (in percent) containing one or 
two photons with a transverse momentum pt{i) > p™"(7) at the Tevatron (pp colhsions at 
y/s = 1.8 TeV). Fractions are obtained with respect to the lowest order cross section. The cuts 
imposed are specified in the text. The relative statistical error on the event fractions from the 
Monte Carlo integration is approximately 1%. 

p??*"(7) (GeV) Z -^ e+e-7 Z -^ e+e--f-f Z -^ e+e-77 [Eq. (||)] 

0.1 52.3 14.6 13.7 

0.3 39.1 7.80 7.64 

1 25.0 3.05 3.13 

3 12.9 0.61 0.83 

10 2.17 1.3 X 10-2 2.4 X 10"^ 

p^^"(7) (GeV) Z -^ p^fi--f Z -^ )U+//-77 Z -^ /i+^-77 [Eq. (|l|)] 

0.1 31.2 4.76 4.87 

0.3 23.8 2.67 2.83 

1 15.4 1.07 1.19 

3 8.19 0.25 0.34 

10 1.62 8.4 X 10-3 1.3 X 10-2 
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FIGURES 




Diagrams by MadGraph 

1 . .41 





graph 1 



graph 2 



graph 3 






graph 4 



graph 5 



graph 6 



FIG. 1. The Feynman diagrams contributing to t'^f^ -^ e^e 77 at tree level, as generated by 
MadGraph. A (Z) represents a photon (Z boson), e an electron or positron, and nujmu a muon 



neutrino, v^. 
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FIG. 2. The distribution of the separation between the two photons, AR^^, for a) pp -^ i~^i~^'j 
and b) pp -^ i'^v^^ at y/s = 1.8 TeV. The sohd and dashed histograms show the differential cross 
sections for electrons and muons, respectively. The cuts imposed are described in the text. 



27 






b 



b 

T3 



1.5 



1.0 



0.08 



a) pp ^1"^1 77 

Vs = 1.8 TeV 

Pt(7)>0.1 GeV 

M 



f 




e 
^ 



"b 
T3 






b 



60 80 100 120 

m(ll) (GeV) 



140 



0.06 



0.04 



0.02 



0.00 



b) pp^l"^i/77 
Vs = 1.8 TeV 
Pt(7)>0.1 GeV 




60 80 100 120 

ititCIp^t) (GeV) 



140 



FIG. 3. Ratio of a) the pp -^ i'^i^^^ and the lowest order pp -^ i'^i^ cross sections as a 
function of the i^i^ invariant mass, and b) the pp -^ i^u^^ and the lowest order pp ^ i'^u cross 
section versus mT{i^x) ^^ V^ — l-^ TeV. The solid and dashed histograms show the cross section 
ratios for electrons and muons, respectively. The cuts imposed are described in the text. 
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